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In  trodur.  Li  on 


Slender  body  theory  has  enjoyed  great  success  in  aerodynamics , and 
in  recent  years  there  has  been  considerable  interest  in  applications  in 
naval  hydrodynamics.  Ship  hulls  are  long  and  slender, and  this  is  a 
logical  feature  to  exploit  in  simplification  of  the  governing  equations. 
Michell's  classical  theory  of  wave-making  resistance  published  in  1898 
considered  the  ship  to  be  a "thin"  body  and  represented  the  hull  by  use  a 
of  distribution  of  "Havelock  sources"  over  a vertical  mid-plane.  How- 
ever, ships  tend  to  be  "slender"  as  opposed  to  "thin"  because  the  beam 
and  draft  tend  to  be  of  the  same  order  of  magnitude  and  both  tend  to 
be  small  in  comparison  to  the  length.  It  was  hoped,  therefore,  that 
slender  body  theory  might  represent  a significantly  better  alternative 
to  Michell's  theory  in  the  area  of  wave  resistance  and  at  the  same  time 
yield  useful  results  for  the  hydrodynamic  forces  in  the  transverse 
direction  both  in  steady  and  unsteady  motion.  Newman  (1970)  has  given 
a review  of  the  application  of  slender  body  theory  to  ship  hulls  which 
discusses  various  kinds  of  hull  motion  and  wave  resistance. 

Major  advances  in  the  application  of  slender  body  theory  in  ship 
hydrodynamics  have  been  made  largely  in  the  past  decade  or  so 
principally  by  Tuck  (1963,1964);  Newman  (1964,1965),  Ogilvie  (1969), 

Orsell ( 1968)  and  Faltinsen(1971 ) . The  method  of  approach  involves 
the  use  of  matched  asymptotic  expansions  which  is  a powerful  method 
for  dealing  with  singular  perturbation  problems.  Potential  flow 
past  a slender  body  is  recognized  as  a singular  perturbation  problem 
and  in  such  cases  an  inner  and  outer  expansion  is  developed.  The 
inner  expansion  is  valid  in  the  near-field,  i.e.,  near  the  body  and 
the  outer  expansion  is  valid  in  the  far-field.  The  unknowns  re- 
maining in  either  of  the  expansion  are  determined  by  a matching 
procedure.  The  general  method  of  matched  asymptotic  expansions 
is  discussed  by  Van  Dyke(1964). 


(1) 


In  the  present  report  the  application  of  the  slender  body  theory 
to  a ship  hull  in  steady  motion  at  high-Froude  number  is  considered. 

While  most  displacement  ships  operate  in  the  zero  to  intermediate  Froude 
number  range, there  is  a new  generation  of  high  speed  ships  generally  of 
the  catamaran  configuration  which  utilize  rather  thin  or  slender  hulls 
and  operate  at  rather  high  Froude  number.  The  surface  effects  ship  hulls 
represent  the  principal  example  of  this  type  of  hull  even  though  the  flow 
about  an  SES  hull  is  somewhat  more  complex  than  the  one  considered  here 
because  of  the  pressure  difference  between  the  two  sides  of  the  hull. 

In  the  following  a single  slender  hull  operating  a high-Froude  number 
is  considered  from  the  view  point  of  slender  body  theory. 

# 

2.  Boundary-value  Problem  For  Ship  Hydrodynamics 


We  assume  the  fluid  to  be  ideal  (inviscid,  incompressible  and 
homogeneous)  and  the  motion  to  be  irrotational . If  (x,y,z)  represents 
Cartesian  coordinates  fixed  in  space  with  the  gravitational  force 
acting  in  the  negative  y-direction  and  t denotes  time,  a velocity 
potential  $(x,y,z,t)  can  be  introduced  such  that  the  fluid  velocity 
vector  is  given  by 


The  condition  of  incompressible  flow  applied  to  (1),  i.e.,7.q  = 0 
yields  Laplace's  equation  as  the  governing  differential  equation, 


V2®  » A + + d>  *0 

xx  yy  zz 

The  fluid  pressure  P(x,y,z,t)  is  expressed  in  terms  of  $ and  the 
elevation  yf  as 


(2) 


P . — p<p t—  *spV<t>.7(J>-pgy  + C(t) 


(3) 


(2) 


where  C(t)  may  be,  in  general,  a function  of  time,  p denotes  the  mass 
density  of  the  fluid,  and  g,  the  acceleration  of  gravity.  The  function 
C(t)  is  at  most  a function  of  time,  and  can  be  determined  by  evaluation 

A A A 

of  P+p^^+^pVft'ViJj+pgy  at  some  reference  point.  This  point  is  generally 
taken  at  the  free  surface  at  infinity  where  all  of  the  values  are  known. 


The  surface  of  the  hull  is  described  by  the  equation  s(x,y,z,t)=0. 
Since  this  surface  is  considered  to  be  impermeable  the  kinematic  boun- 
dary condition,  which  imposes  the  condition  that  the  relative  fluid  velo- 
city in  a direction  normal  to  the  hull  must  vanish,  is  given  by 


Ds 

Dt 


• Vs-0 


on  s-0 


(4) 


On  the  free  surface  there  are  two  boundary  conditions  which  must 
be  applied.  If  the  elevation  of  the  free  surface  is  described  by 
y«n(x,z,t)  the  kinematic  boundary  condition  obtained  in  a fashion 
similar  to  (4)  results  from  setting  the  substantial  deriative  of  the 
surface  n-y-0  to  zero, 


-^(n-y)- 


■n  +$  n +<t>  n -$ 
t x x z z y 


0 on  y-n 


(5) 


In  addition  to  (5)  a dynamic  boundary  condition  must  also  be  imposed. 
This  boundary  condition  results  from  simply  applying  P-0  on  the  free 
surface  and  is  given  by 

A A A 

*V4i+gy-C(t)  on  y-n  (6) 

where  C(t)  is  evaluated  with  the  aid  of  a known  reference  point  in  the 
fluid. 


(3) 


In  addition  to  the  boundary  conditions  (4-6)  the  potential  function 
must  satisfy  certain  radiation  conditions  at  x2+y2  -*•  00  appropriate  to 
the  particular  problem  under  consideration.  Also,  on  the  bottom  the 
kinematic  boundary  condition 


0 on  y»- h 


(7) 


is  applied  if  the  water  depth  is  finite  and  denoted  by  h. 


3.  Perturbation  Procedure  For  a Slender  Body 

At  this  juncture  let  us  consider  a particular  problem:  a slender 
hull  is  positioned  with  its  long  axis  coincident  with  the  x-axis  and 
held  fixed  in  space  as  indicated  in  Figure  1.  A uniform  stream  of 


Figure  1 Definitions 

magnitude  U makes  an  angle  a with  respect,  to  the  x-axis  so  that  the 
resulting  x,y,  and  z components  are  -U  cos  a,  0,  -U  sin  a,  respectively. 
The  coordinate  system  is  placed  at  the  mean  free  surface  so  that  the 
y«0  plane  represents  the  mean  free  surface, and  the  water  depth  is  taken 
as  infinite.  The  motion  will  be  restricted  to  steady  motion, and  the 


(4) 


hull  Is  described  by  s(x,y,z;e)  * 0 where  the  parameter  e * beam  or 
draft -to -length  ratio. 

Outer  Problem 


The  problem  which  is  here  defined  is  recognized  to  be  dependent 
on  two  dimensionless  parameters,  the  beam/ length  ratio,  e,  and  the 
incidence  angle,  a.  Thus,  as  an  initial  step  we  may  expect  the  solu- 
tion to  be  expandable  in  a power  series  in  both  c and  a.  (At  least 
the  solution  should  behave  in  a regular  manner  with  the  exception  of 
the  region  near  the  body.)  Introducing  the  body  potential,  $,  the 
total  potential  may  be  expressed  as 

$ * <t>(.x,y,z;e,a)-Ux  cos  a-Uz  sin  a (8) 

where  the  body  potential  4>  is  expressed  by 

$=e$j+a4>“+O(e2)+0(a2)  (9) 

in  the  outer  region.  The  first  term  of  the  Taylor  series,  (0,0)=  <t>0,  is, 
of  course,  zero  because  if  both  e=0  and  a=0  the  uniform  stream  is 
directed  parallel  to  a line  of  zero  thickness  and  the  body  potential 
vanishes.  Thus,  expanding  sina  and  cosa  about  a=0  in  (8)  and  substi- 
tuting into  (9)  gives  the  total  potential  as 

$ ■ e<^  + a($^-Uz)-Ux+0(e2)-K) (a2)  (10) 

and  substituting  this  form  of  the  potential  into  (2)  shows  that 

72*j»V2*“  =0  (11) 

since  a and  e are  independent  parameters. 


(5) 


In  similar  fashion, we  may  express  the  free  surface  elevation  in 
the  form  of  (9), 

n“en^an^+°(  e2)+0(a2)-f  O(ac)  (12) 

Substituting  (8) , (9)  and  (12)  into  (5)  and  collecting  coefficients  of 
like  powers  of  e and  a gives: 

<f>?  + Unf  =0  on  y=0  (13) 

ly  lx 

and 

+ Un?  = 0 on  y=0  (14) 

ly  lx 

The  dynamic  free  surface  boundary  condition  is  obtained  by  use  of 
(6)  by  evaluating  C(t)*liu2  at  x=°°  on  y*0  and  substituting  in  the  form 
of  $ and  n as  given  in  (12)  and  (8).  Here  again  two  separate  boundary 
conditions  result  because  of  the  independence  of  the  two  small  parameters, 
c and  a.  The  first  of  these  is  obtained  by  setting  a=0  as 

nf+  — 4>f  =0  on  y*>0  (15) 

1 g rlx 

and  the  second  by  setting  £”0, 

n?+_ * 0 on  y*0  (16) 

1 g lx 


f Eqs.  (13)  and  (14)  together  with  (15)  and  (16)  yield  the  linear  free 

surface  boundary  conditions. 


Eqs.  (17)  and  (18)  represent  familiar  linear  boundary  conditions  for 
the  free  surface.  If  we  now  let  the  Froude  number  U2/gL  (where  L Is 
considered  to  be  0(1))  become  large,  (17)  and  (18)  become  simply 


on  y=0 


(19) 


“ 0 


on  y*0 


(20) 


In  summary,  we  have  derived  a boundary-value  problem  for  the  outer 
potential.  The  solution  is  singular  in  the  region  of  the  body  so  the 
hull  kinematic  boundary  condition  is  not  applied  to  the  outer  solution. 
We,  therefore,  have  in  summary  for  <f>j. 


*lyy+  *lzz’  ° 


* 0 on  y«0 


(22) 

(23) 


and  similarly  for  $ , 


a a a 
1 + *,  + *7 


•XX 


lyy 


izz 


(24) 


* 0 on  y»0 


(25) 


The  boundary  condition  on  the  hull  is  in  this  case  to  be  replaced  by  a 
matching  condition. 


Inner  Problem 


We  now  consider  the  inner  potential  which  is  valid  near  the  slender 
body.  The  solution  to  the  slender  body  problem  is  singular  in  the  inner 
region  but  the  degree  of  singularness  at  the  outset  is  unknown.  We  there- 


(7) 


fore  begin  by  writing  the  body  potential  as 


$ - A0 (e)^ (x.Y.ZjaJ+i! (x,Y,Z, ;a)+. . . (26) 

where  the  Y and  Z coordinates  are  stretched  by  the  parameter  e, 

Y=y/e  (27) 

Z=z/e  (28) 

and  Ao(e),A1(e)>  etc.  are  called  gage  functions  as  defined  by  Van  Dyke 

(1964).  These  gage  functions  will  eventually  be  determined  by  use  of 
the  principle  of  least  degeneracy  also  discussed  by  Van  Dyke. 

The  body  potential  is  expressed  in  (26)  as  an  inner  perturbation 
series  in  e,  the  beam-length  ratio.  Here  the  incidence  angle,  a,  is 
allowed  to  remain  finite  and  is  not  considered  in  any  limiting  process. 
We  shall  return  to  this  point  later. 


If  the  inner  representation  of  the  body  potential  given  by  (26) 
is  substituted  into  (8)  and  the  result  substituted  into  the  Laplacian, 
(2),  we  find. 


+<t>  +4>  + Al^ 

°xx  °YY  °ZZ  A2(e) 


Te2^  + 

L kx 


+ t*  i 


YY 


1 1 + 

ZZJ 


. ..  = 0 (29) 


Clearly,  if  A^eJ/Agfc)  is  at  least  of  0(e)  then  the  three  dimen- 
sional Laplacian  degenerates  to  the  two-dimensional  form  for  the  leading 
term  in  the  inner  expansion: 


$ + $ »o  (30) 

°*Y  °ZZ 

The  inner  representation, being  valid  near  the  hull .must  satisfy  the 
boundary  condition  on  the  hull.  If  s is  expressed  in  inner  variables  as 


(8) 


s(x,y,z;e)  = S(x,Y,Z;e)  = 0 then  the  hull  kinematic  boundary  condition, 
which  is  obtained  uy  use  of  (26)  and  (8)  with  (4). gives 


Sx+* 


oyV*° 


Al  (g) 

A (c) 
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i S +* 

X X 


1 SY+Ii‘ 

Y Y 


e2USj{  _ eUaSz 

A (e)  A (e) 

o o 


(31) 


Now  if  a were  to  be  considered  to  be  fixed  at  0(1)  then  the  principle  of 
least  degeneracy  would  require  that  Aq(c)  = e,  and  the  hull  kinematic 
condition  would  read 


VYHoZSZ*,faSZ  °"  S-° 


(32) 


This  boundary  condition  simply  indicates  that  in  any  transverse  plane 
perpendicular  to  the  x-axis  the  body  potential  must  just  cancel  the 
cross-flow  velocity,  Ua.  It  is  noted  that  in  such  a case  the  cross- 
flow  is  an  order  of  magnitude  larger  than  the  perturbation  associated 
with  the  thickness  distribution  of  the  body  and  thickness  effect  is, 
therefore,  lost. 


Alternately,  if  the  angle  of  attack  were  taken  to  be  zero  then  the 
principle  of  least  degeneracy  leads  one  to  A^fe)-?:2  and  (31)  gives. 


4>  S +$  S -US  on  S-0 
OY  Y OZ  Z X 


(33) 


a familiar  result  in  slender  body  theory. 


A third  possibility  is  to  take  a to  be  of  the  same  order  of  magni- 


(9) 


tude  as  e and  retain  both  of  the  terms  on  the  right  hand  side  of  (32) 


and  (33).  This  gives  Aq(e )«e2  and 

S V*o.;SZ  ' U(Sx  + f V °n  S“° 
An  alternate  way  of  writing  (34)  is 

- &*0"m  a \ 

TN"  u\Sn  e 3Z/ 


(34) 


(35) 


where  N»n/c  and  n represents  distance  measured  normal  to  the  curve  in  the 
Y-Z  plane  described  by  the  equation  S(x,Y,Z;e)  = 0 with  X held  fixed. 


The  free  surface  condition  to  be  applied  to  the  inner  solution  is 
obtained  by  eliminating  n between  (5)  and  (6).  The  boundary  condition 
obtained  in  this  way  is 

$^(7$ • 7$^)+^ (^<f> ’^$z)+g'?y=0  on  Y**n/e  (36) 
Substitution  of  (26)  and  (8)  into  (36)  gives 


(e2$0x-u)  [e2*0xx^^-U)+e2*oY*oxy+e«0xZ  <£*o 

+(e*oZ-Ua)  [e*oJe\x-U)+c%Y%YZ+%ZZ(£%Z-Ua)+g£\m°  °n  Y’° 

Now,  if  e-»0  with  a-fixed,  (37)  gives 


(37) 


*0zz”°  °n  Y"°  (38) 

or  from  this  we  can  deduce  that 

4>  -0  on  Y-0  (39) 

o 

is  the  appropriate  boundary  condition  in  the  inner  region.  The  problem 
for  the  first  term  of  the  inner  expansion  for  the  high-Froude  number 


(10) 


slender  body  theory  is,  in  summary 


<t  +♦  “0 

oYY  oZZ 


3»o  „,sx  jO  3N  . _ 

ti?  *7  az  ) on  s=0 


(40a) 


(40b) 


<t>0  =0  on  Y«0  (40c) 

There  is  no  radiation  condition  specified  in  (40);  this  condition  is 
replaced  by  matching  with  the  outer  solution. 


4.  Solution  for  the  Inner  Problem 

The  inner  problem  specified  in  (40)  is  essentially  a two-dimensional 
problem  as  indicated  in  Figure  2.  Such  problems  can  be  solved  for  rather 


Figure  2 Inner  Problem 

arbitrary  shapes  by  use  of  a Green's  function  of  the  form 


G(Z,Y;c,n)-ln^(Z-02+(Y-n)2-lnV(z“^2  +<Y+h>2 


(ID 


The  first  term  in  (4 1) represents  a source,  and  the  second  term  represents 
a negative  image  with  respect  to  the  Z axis.  The  potential  ^may  be 
represented  by  a surface  distribution  over  the  hull  of  these  two-dimen- 
sional singularities. 

In  particular,  the  potential  may  be  represented  by 


o 2 it  J F(x;n,OG(Z,Y;e,n)dC 


where  C denotes  the  curve  in  the  Y-Z  plane  representing  the  cross- 
section  of  the  hull  and  dC  denotes  an  arc  length  along  that  curve. 
Using  (41) for  G,  (40a)  and  (40c)  are  automatically  satisfied.  The 
hull  boundary  condition  (40b)  then  yields  the  integral  equation  for 
the  source  strength  distribution,  F(x;Y,Z)  as 


F(x;Y,Z)+^  /F(x;n,0  (Z,Y;c,n)dC 
rr  u N 

-Km)  °”c 


This  equation  can  be  solved  numerically  by  subdiving  C into  a finite 
number  (say  N)  of  segments.  Eq.  (43)  can  then  be  reduced  to  a set  of 
N linear  equations  in  N unknown  values  of  F^  over  each  of  the  segments. 


Cylindrical  Hull 


For  the  case  of  a semi-immersed  circular  cylinder, a multipole 
series  may  be  used  to  represent  the  solution.  Consider  the  surface 
described  in  cylindrical  coordinates  by 


S(R,0,x)-R-RH(x) 


where  R-r/e  denotes  the  radial  coordinate  in  inner  variables  and  R^ 


(12) 


denotes  the  radius  of  the  hull  in  the  inner  variable  (R^  -r^/e)  and 
is,  in  general,  a function  of  x.  The  kinematic  boundary  condition  on 
the  hull,  (40b),  takes  the  form, 

||°  - UO^  + 2.  cos0)  on  R-RR  (45) 

where  the  polar  coordinate  system  is  related  to  the  cartesian  system  by 


Z-R  cos  0 
Y— R sin  0 


(46a) 

(46b) 


The  solution  for  $ satisfying  both  (40a)  and  (40b) is  given  by 
o 


VE  an  3in(n0) 


n-1 


Application  of  the  hull  boundary  condition  (45)  now  gives 
oo  a n 

E rfen  sin  n9-URH  - f u cos0 

n-l  RH 


(47) 


(48) 


By  use  of  the  results 
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we  find 
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or,  carrying  out  the  integration. 
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The  inner  potential  function  is,  therefore,  represented  by  the  series 


s-  £ (=fcr)  S' 
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(56) 
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RHn+1  sin(n9) 


The  dynamic  pressure  on  the  hull  can  be  obtained  from  Bernoulli's 
equation,  disregarding  terms  of  0(e4)  and  o(e2a2), 
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P-0  e2 


U® 


'°x"  (*°r  + R7  *o0)“(^)^R*oesine+(f)  U2sin2ej  (57) 


The  normal  (or  yaw)  force  per  unit  length  along  the  hull  can 
be  determined  through  integration  of  the  pressure  over  the  immersed 
surface.  That  is, 

’ it 

P RH  cos9  doj  dx  (58) 


(14) 


After  substitution  of  (56)  into  (58)  and  some  lengthy  manipulations 
and  integrations  we  get,  for  the  normal  force  on  the  segment  dx, 

dV  *c2aPu2  7 RH  rh'  2+2  - n^+(1~1)-  ^ dx  (59) 

The  net  normal  force  is  obtained  by  integration  from  the  bow  where 


Rh=0  to  the  stern  where  , 


F«-  i pu2  f r2HS  h 2+X  r 


-(-l)n  n(l+(-l)n) 


The  dimensionless  normal  (or  yaw)  force  coefficient  based  on  the  stern 


area,  therefore,  becomes 


2 SeU2*  r2 


77  (2+^ 


l-(-l)n  n( l+(-l)^) \ ( _ (61) 

n 


where  r is  defined  as  =r„  /e  and  rH  represents  the  hull  radius  in  outer 

HS  Hs 

variables.  For  boat-tailed  hulls  we  get  the  familiar  result  Fz»0  because 


RHc”0- 


The  summation  in  (61)  can  be  somewhat  rearranged  to  read 


l n-3,5,7  1 


The  series  converges  rapidly  to  approximately  0.128  so 


f -0.75a  (63) 

z 

The  value  of  f2-0.75a  for  the  cylindrical  hull  at  high  Froude  number  can 


be  compared  with  the  classical  value  of  f -2a  for  the  same  hull  operating 


5 


with  the  low-speed  free  surface  boundary  condition,  34>  /3Y=0  on  Y=0. 

o 

It  is  reasonable  to  expect  a somewhat  smaller  value  for  fz  for  the 
high-Froude  number  case  because  the  high-Froude  number  free  surface 
boundary  condition  makes  the  Y=0  plane  act  like  a zero  pressure  sur- 
face and,  therefore,  right  at  the  free  surface  there  can  be  no  pres- 
sure difference  across  the  hull.  The  low  Froude  number (actually  low 
to  moderate)  boundary  condition  makes  the  free  surface  act  like  a rigid 
wall  and  a pressure  difference  across  the  hull  can  exist  at  the  Y-0 
level. 


The  moment  about  the  Y-axis  is  given  by 

b 


x»-a 


(64) 


where  a and  b denote  the  distance  from  the  origin  to  the  stern  and  bow, 
respectively, so  that  the  hull  length  is  Lma+b.  Thus,  using  (59)  with 
(64)  gives 


l-(-D 
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n(l+(-l)n 

n2-l 


(65) 


where 

As-area  of  the  stem  section  below  the  mean  waterline 
V”hull  displacement 
L**hull  length 


When  the  series  is  evaluated  we  get 


(66) 


(16) 


and  for  boat-tailed  (or  closed)  bodies  the  moment  coefficient  becomes 
n»y*0.75a.  This  value  may  be  contrasted  with  the  low  Froude  number 
result  of  2.0a  which  is  often  referred  to  as  Munk's  moment. 

5 . Conclusions 

Slender  body  theory  has  been  applied  to  a slender  ship  hull  opera- 
ting at  high-Froude  number.  For  the  case  of  a semi-circular  hull, closed- 
form  results  were  obtained  for  the  yaw  force  and  moment.  These  values 
are  3 ess  than  half  the  values  of  the  equivalent  coefficients  associated 
with  the  rigid-wall  free  surface  boundary  condition. 


(17) 
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